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Abstract 

Based on the phenomenological Skyrme interaction various density-dependent nuclear 
matter quantities are calculated up to second order in many-body perturbation theory. The 
spin-orbit term as well as two tensor terms contribute at second order to the energy per 
particle. The simultaneous calculation of the isotropic Fermi-liquid parameters provides a 
rigorous check through the validity of the Landau relations. It is found that published results 
for these second order contributions are incorrect in most cases. In particular, interference 
terms between s-wave and p-wave components of the interaction can contribute only to 
(isospin or spin) asymmetry energies. Even with nine adjustable parameters, one does not 
obtain a good description of the empirical nuclear matter saturation curve in the low density 
region 0 < p < 2po. The reason for this feature is the too strong density-dependence of 
several second-order contributions. The inclusion of the density-dependent term is 

therefore indispensable for a realistic description of nuclear matter in the Skyrme framework. 


1 Introduction and summary 

Mean-field theories are widely used to analyze a large variety of many-body systems for which the 
independent particle picture can be applied. In nuclear physics, mean-held approaches generally 
lead to satisfactory results for bulk properties of atomic nuclei, such as their masses, radii and 
ground-state deformations. However, in some cases improvements beyond the mean-held treatment 
are dehnitely needed. For instance, mean-held models do not predict accurately the single-particle 
spectra (i.e. energies, spectroscopic factors and fragmentations of single-particle states), especially 
close to the Fermi energy. A possibility for improvement is given by introducing the coupling 
between individual nucleon degrees of freedom and collective degrees of freedom. 

Phenomenological Skyrme forces have been used widely and successfully for non-relativistic 
nuclear structure calculations of medium-mass and heavy nuclei niEi. In ehect such purely phe¬ 
nomenological two-body interactions can be viewed as to provide a convenient parameterization 
of the nuclear energy density functional on which the self-consistent mean-held treatment (in the 
form of solving the Kohn-Sham equations) can be based. In order to go beyond the common 
mean-held (Hartree-Fock) approximation, second-order calculations with the phenomenological 
Skyrme force have been pursued in recent works by Moghrabi, Grasso et al. P in El E]. These 
exploratory studies have been concerned with the second-order contributions to the equation of 
state of inhnite isospin-symmetric nuclear matter and pure neutron matter. Since the (two-body) 
Skyrme interaction is of zero range, divergences occur at second order in many-body perturbation 
theory, which have been treated either by cutoh regularization Pi or by dimensional regulariza¬ 
tion p. It has been concluded in refs. iP that after readjusting the parameters of the Skyrme 
interaction a simultaneous good reproduction of the equation of state of isospin-symmetric and 
isospin-asymmetric nuclear matter can be achieved at second order. 

However, a closer inspection of the published expressions for the second-order contributions in 
ref. P raises doubts about their completeness and correctness. The hrst dehcit is the omission of 
a second-order contribution from the spin-orbit term proportional to Wq. This p-wave interaction 
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term averages to zero in infinite matter at first order, but not at second order (see e.g. appendix 
A in ref. [7] for the second-order contributions from a “boson-exchange” spin-orbit interaction). 
Secondly, the stated expressions |5] for the energy per particle of isospin-symmetric nuclear matter 
and pure neutron matter include interference terms between s-wave and p-wave components of the 
Skyrme interaction (i.e. terms proportional to and ^ 1 ^ 2 )- If is very questionable how such a 
mixing of different parities could be induced by the (homogeneous) nuclear medium. Moreover, 
the expressions for the cutoff-independent second-order contributions given in a later paper 0 by 
the same authors disagree partly with those published in ref. [5]. 

These apparent inconsistencies provide good reason to reconsider and extend the calculation 
of the Skyrme interaction to second order in nuclear matter. The present paper is organized as 
follows: In section 2, the Skyrme force is introduced as a Galilei-invariant two-nucleon contact- 
interaction in momentum space. On this basis, hrst-order Hartree-Fock results are rederived for 
a variety of density-dependent nuclear matter quantities, such as the energy per particle and the 
(isospin and spin) asymmetry energy. Section 3 is devoted to the second-order calculation in many- 
body perturbation theory. We employ a scheme where the in-medium (particle-hole) propagator is 
decomposed into the vacuum propagator and a medium-insertion. This scheme has the advantage 
that ultraviolet divergences appear only in the second-order Hartree and Fock diagrams with 
two medium-insertions. The pertinent one-loop rescattering in vacuum is treated by dimensional 
regularization and shown to vanish identically. As a consequence the genuine second-order terms 
arise entirely from the Hartree and Fock diagrams with three medium-insertions. The occurring 
finite integrals over three Fermi spheres can be solved analytically and the results for the nuclear 
matter quantities consist of even powers of the Fermi momentum kf multiplied with products of 
the Skyrme parameters and a numerical coefficient of the form ln2 -|- r^, where Vj is some rational 
number. It is found that the energy per particle of isospin-symmetric nuclear matter E{kf) and 
pure neutron matter En(kn) do not include any s-wave and p-wave interference terms at second 
order. However, such sp-interference terms show up in the (isospin and spin) asymmetry energies. 
The reason and condition for sp-mixing to emerge is explained in appendix A by considering the 
in-medium loop for a proton-neutron pair in nuclear matter with unequal proton and neutron 
densities. In section 4, the effective nucleon mass and the isotropic Fermi-liquid parameters are 
calculated up to second order. These results allow for a rigorous check through the validity of 
several Landau relations [S]. In section 5, the effects of tensor interactions at second order are 
computed. Actually, with the inclusion of two independent tensor terms the Skyrme force gets 
completed to the most general two-nucleon contact-interaction quadratic in momenta. In section 
6, the iteration of the Skyrme force with the one-photon exchange is considered. This mechanism 
generates a charge-symmetry breaking contribution to the energy density of protons, which could 
help to resolve the Nolen-Schiffer anomaly [5]. This anomaly refers to the observation that the 
binding energy differences of mirror nuclei cannot be explained by the Coulomb interaction of the 
protons alone. Remarkably, one finds an infrared-finite result for the energy per proton Ep{kp) 
even though a regulator photon mass has be introduced initially. Section 7 deals with some 
provisional parameter fits. It is found that with nine adjustable parameters in the complete 
second-order expression, one does not obtain a good reproduction of a realistic nuclear matter 
saturation curve in the low-density region 0 < p < 2po- The reason for this failure lies in the 
too strong density-dependence of several second-order contributions. From this observation 
one can conclude that the inclusion of (ad hoc) density-dependent terms is indispensable for 
a realistic description of nuclear matter in the Skyrme approach. At a formal level, these are 
taken into account by performing appropriate parameter substitutions, to to + and 

toXo —)• toXo + Itsa^sp^/®, in all first- and second-order expressions derived in this work. 

Let us clarify that the main purpose of this paper is to compile the (correct) analytical ex¬ 
pressions for various nuclear matter quantities as they arise from a second-order calculation with 
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a two-body contact-interaction of the Skyrme type. The question of whether improvements are 
possible in comparison to a hrst-order mean-field approach has to be explored in future studies. 

2 First-order Hartree-Fock calculation 

This section serves mostly for orientation. We introduce various density-dependent nuclear matter 
quantities and specify their hrst-order Hartree-Fock contributions. The Skyrme interaction is a 
zero-range two-nucleon contact-interaction with transition matrix-element in momentum-space: 

14k = 4(1 ~\-XQPa) 4^(1 + *?in) + 4(1 + 3:2T’cr) 5'out ‘ 5'in + * hFo(cTi + (72) ' (<?out X Qin) ■ (1) 

Here, Po- = (1 + <?i • <?2)/2 is the spin-exchange operator and gin = (pi —p2)/2 and ^ut = {pi —pi) 12 
denote half of the momentum difference in the initial and hnal state, respectively. This particu¬ 
lar form of Vsk, involving only external momentum differences, is dictated by Galilei invariance 
(and invariance under time-reversal: ^ —^jut, <?i ,2 —t —^ 1 , 2 )- The individual terms in eq.(l) 

stand for: s-wave contact interaction (to), s-wave effective range correction (ti), p-wave contact 
interaction (^ 2 ), and spin-orbit interaction (Ho). The general contact-potential up to order p^, as 
it is commonly used in chiral effective held theory UniElDE], has two additional tensor terms 
which will be studied separately in section 5. Note that we have not explicitly written in eq.(l) 
the density-dependent Skyrme term 14(1 + Formally, it can be included through the 

substitution of parameters 4 to + and fo^^o IqXq -|- all following formulas. 

However, from a rigorous point of view such density-dependent couplings will lead to a conhict 
with basic many-body relations (Landau relations [8] and Hugenholtz-Van-Howe theorem |13j). 
whenever derivatives with respect to density p are involved. 



Figure 1: Upper row: First-order Hartree and Fock diagrams with two medium-insertions and 
combinatoric factor 1/2. The wavy line represents the Skyrme interaction Vsk. Lower row: Second- 
order Hartree and Fock diagrams with three medium-insertions and combinatoric factor 1. 

The hrst-order Hartree and Fock diagrams generated by Vsk are shown in the upper part of 
Fig. 1. The short double-line symbolizes the medium-insertion (i.e. the diherence between the in¬ 
medium and vacuum nucleon-propagator) proportional to the step-function 6{kf — \pj\). The Fermi 
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momentum kf is related the (total) nucleon density in the usual way, p = 2kj/37i‘^. Performing 
the appropriate spin- and isospin-traces and integrals over two Fermi spheres, one obtains for the 
energy per particle of isospin-symmetric nuclear matter the well-known hrst-order expression: 

E(k,y^“^ = + h [3(i + *2(5 + 412 )] I, (2) 

In the same way, but leaving out the isospin-trace, the energy per particle of pure neutron matter 
comes out as: 

= ^{|(1 - + i b(l - + 3i'2(l + ^ 2 )] } . (3) 

with kn the neutron Fermi momentum related to the neutron density by pn = k^/37i‘^. Another 

important quantity of high interest is the isospin asymmetry energy of nuclear matter. An isospin- 
asymmetric conhguration with unequal proton and neutron Fermi momenta kp^n = /c/(l =F is 
realized by the substitution: 

- \pj\) ^-^e{kp - \pj\) + ^-^e{kn - \Pj\) , (4) 

for the medium-insertion and an additional expansion of the (asymmetric) energy per particle: 

EUkp, kn) = E{kf) + S^A{kf) + 0(6^) , (5) 

around 5 = 0 dehnes the density-dependent isospin asymmetry A{kf) of nuclear matter. The 
hrst-order Hartree-Fock contribution from the Skyrme interaction Vsk has the well-known form: 

[^ 2(4 + 5 x 2 ) - Stixi] I. (6) 

In an analogous way the density-dependent spin asymmetry energy S{kf) of nuclear matter can 
be constructed. The substitution in spin-space: 

^{kf - \Pj\) - \Pj\) + ^-^G{ki - \Pj \), ( 7 ) 

with k-f^i = kf{l ± 77 )^/^ the Fermi momenta of the spin-up and spin-down nucleons, is followed by 
an expansion of the energy per particle of spin-polarized nuclear matter around 77 = 0 : 

4oi(A:t, %) = E{kf) + p^S{kf) + 0{p ^). (8) 

In the absence of any quantitative empirical knowledge, the minimal condition to be satished by 
the spin asymmetry energy is that it must be positive: S{kf) > 0. This guarantees the spin- 
stability of nuclear matter by excluding a phase transition to a ferromagnetic state of lower energy 
density. The hrst-order Hartree-Fock contribution to S{kf) from the Skyrme interaction Vsk reads: 

= ^ 2^2 |^o(2a;o ~ [StiXi + ^2(4 -|- 5 x 2 )] j- . (9) 

One notices that the expression for can be obtained from A{kf)^^^^'^ in eq.( 6 ) by just 

changing the signs of the two parameters Xq and Xi. The reason behind this connection is an equi¬ 
valent form of the contact-interaction which results from Vsk in eq.(l) through the substitutions: 

to(l + a^o-Po-) to{l — XqPt) , 
tl{l+XlPa) tl{l — XlPr) , 

0 (l+a^ 2 -P(T) f 2 (l + a^ 2 -Pr) ; ( 10 ) 
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with Pr = {1 + ti-T2)/2 the isospin-exchange operator. The known eigenvalues of and Pr are 
(—and (—1)'^^^ with S,I E {0,1} the total spin and isospin quantum numbers. Since the 
parameters to and ti act only in s-wave states with S + I = 1, one has effectively for s-waves 
the equality P„ = —Pr- Likewise, the parameter t 2 acts only in p-wave states with S' -|- / = 0, 2 
and therefore for p-waves the opposite equality Po- = Pr holds. Another obvious consequence 
of the equivalent interaction is that only the combinations to(l — a;o), ti(l — xi) and t 2 (l + X 2 ) 
can contribute in neutron matter, since P^ = 1 for a two-neutron state. This is exemplified by 
calculating the first-order contribution to the spin asymmetry energy of pure neutron matter: 

= ^2^2 ~ ~ 1) + 3^2(1 + 2 ^ 2 )] j" • (H) 


Again, the density-dependent quantity Sn{kn) must be positive: Sn{kn) > 0, in order to prevent 
neutron matter from a ferromagnetic instability at low and moderate densities. We conclude this 
section by giving the contributions from the relativistically improved kinetic energy: 




lOM 


A{kfY^'^'> = s{kfY^'^'> = 


where M = 939 MeV denotes the average nucleon mass. 


3kj 

EniknY^'^^ = 

3kl 

3ki 

(12) 

56M3 ’ 

' lOM 

56M3 ’ 


Sn ( kn, 

j(kin) _ 

kl 

kt 

(13) 

6M 12M3 ’ 


QM 

12M3 ’ 


3 Second-order Hartree-Fock calculation 


In this section we derive analytical expressions for the contributions to the introduced nuclear 
matter quantities as they arise from a second-order calculation with the Skyrme interaction Vsk- 
At second order in many-body perturbation theory the energy density is represented by closed 
three-loop Hartree and Fock diagrams with either two or three medium-insertions. The first set 
involves the rescattering of two nucleons in the vacuum (i.e. without Pauli blocking corrections). 
For the contact-interaction Vsk this one-loop diagram leads to principal-value integrals of the form: 


dH M 
(27r)3 12 _ ^ 



(14) 


with q = (pi — P2)/2, where pi ^2 are two external momenta inside the Fermi sphere |pi^ 2 | < kf. 
The first factor in eq.(14) is the energy denominator and the other factors in curly brackets 
come from the Galilei-invariant contact-interaction. After angular integration and partial fraction 
decomposition of the radial integrand one gets power divergences of the form fY°dl{l, l^, l^}, which 
are set to zero in dimensional regularization. The remaining convergent principal-value integral 
dlij? — q^)~^ = 0 vanishes identically. As a side remark we note that cutoff regularization would 
provide additional terms A^’^’^, and A(g^)^ proportional to odd powers of the cutoff A. For 

the first two the A-dependence could be absorbed into the coupling constants tj and tjXj, while 
the last one would require the introduction of a four-derivative contact-interaction. 

In dimensional regularization (or after renormalization) one is therefore left with the Hartree 
and Fock diagrams with three medium-insertions shown in the lower part of Fig. 1. These lead to 
(finite) principal-value integrals over three Fermi spheres of the form: 


d?pid?p2d?pz M 


\pj\<kf 


(2vr) 


qi-q2 


H{qi,qi-q2,ql. 


(15) 
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with qi ^2 = P 3 — Pi, 2 - The hrst factor M/{qi-q 2 ) is again the energy denominator and the function 
is a polynomial of total degree two or less in its three variables. In order to evaluate eq.(15), 
one eliminates p 2 in favor of ^ and integrates over a shifted Fermi sphere \q 2 — Psl < kf. The 

upper boundary for the dq' 2 -integration is then p^y + ~ P 3 + with the directional cosine 

y = Q 2 -P 3 - By applying this procedure to evaluate the nine-dimensional principal-value integrals, 
one is able to obtain analytical results for all second-order contributions. Let us also note that for 
s-wave interactions ~ to, ti the Hartree and Fock diagram contribute with opposite sign, whereas 
for p-wave interactions ~ t 2 , Wq the Hartree and Fock diagram contribute with equal sign. 

Putting all pieces together one ends up with the following expression for the second-order 
contribution from Vsk to the energy per particle of isospin symmetric nuclear matter: 


E{kf) 


(2nd) _ 


Mk^f f „ 

/ J Q4-2 


2e, 


+ i3(ll - 2 In 2) + + x„ii)(167 - 24 In 2) 

k^ 

^ tl{l + a;^)(4943 - 564In2) + + 8 x 2 + 5a;^)(1033 - 156 In2) 


396 


1188 


99 


- 102 In 2 ) \ . 


(16) 


The hrst important observation one can make is that there are no interference terms between s-wave 
and p-wave components of the Skyrme interaction Vsk. A detailed comparison to the expressions 
written in eqs.(16,17) of ref. [5] reveals that therein only the term proportional to ^^(1 -f Xq) is 
correct. The different form = 3Mkj{Ao + A^kj -f A^k'^^)/2n‘^ given in eq.(5) of the 

later paper [ 6 ] by the same authors represents no improvement either. It be should noted that 
ref. |5] did not use the advantageous scheme of expanding in the number of medium-insertions and 
therefore divergent integrals of the form eq.(14) as well as hnite integrals of the form eq.(15) had 
to be continued together to d 7 ^ 3 space-dimensions. Moreover, the Galilei-invariant formulation 
of the interaction Vsk in momentum-space is important. 

By adding half of the (previous) Hartree contribution to the Fock contribution one arrives at 
the energy per particle of pure neutron matter: 


As expected from the consideration of the equivalent interaction in section 2 , the second-order 
result involves only the parameter combinations fo(l ~ 2 : 0 ), ti(l — Xi) and ^ 2(1 + X 2 ). 

The corresponding expression given in ref. [5] is correct only with respect to the term proportional 
to toti{l — xo)(l — xi), and the (simplest) term proportional — xoY misses a factor 2 therein. 
Note that the factor (1033 — 1561n2)/385 appears also in eq.(32) of ref.[Tl], where the second- 
order contribution to En{kn) has been extracted from the resummation to all orders of particle-hole 
ladder diagrams generated by a spin-independent p-wave contact interaction. 

The second-order contribution to the isospin asymmetry energy A{kf) is more difficult to 
calculate. One has to evaluate now second derivatives with respect to <5 at d = 0 of principal-value 
integrals over three Fermi spheres with different radii kp^n = fc/(lT<5)^^^. The complete expression 
for the second-order contribution from the Skyrme interaction Vsk to the isospin asymmetry energy 


r op2 

= ^ 1 ^ 0(1 - ^o)^(ll - 21n2) + - Xo)(l - a:i)(167 - 241n2) 

+ Y^t?(l - a;i)2(4943 - 564ln2) + ^4(1 + a;2)'(1033 - 156 ln2) 


8fcn 

'297 


W^{<o3l - 102 In2) 


(17) 
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of nuclear matter reads: 


,s Mki f o r o 1 

’ = ^|«o[(l+4)(l-21n2)-10i„] 

+ ^toti [4(1 + a;oa;i)(51 - 221n2) + 7{xo + a;i)(41n2 - 57)] 

k^ r /1301 \ 

+^tl {l + xl)( -461n2) +a;i(481n2 - 607) 

63 L V 6 / 

[ 5(1 + 2 ^ 2 ) -78 In 2 ) +a; 2 (3187-624 In 2) 

k'j 

+ —^ 0 ^ 2(2 Xq x -2 -\- 2 x 03 ^ 2 ) (4 In 2 — 7) 

6 

2k j 

+ -^tit2{2 + xi+ X2 + 2xiX2){l2\n2 - 31) 

63 

4^4 'V 

+Yg5M'o(761 - 1321n2)|. (18) 

Note that the third-last and second-last contribution proportional to to ^2 and tit 2 are sp-mixing 
terms. These arise solely from the second-order Hartree diagram and they are made possible by 
the different radii of Fermi spheres involved. The reasons and conditions for sp-interference to 
emerge for a nuclear matter quantity are explained in appendix A by considering the in-medium 
loop for a proton-neutron pair in nuclear matter with unequal proton and neutron densities. 

In the same way (i.e. by evaluating derivatives 8“^/dTf at p = 0) the second-order contribution 
to the spin asymmetry energy S{kf) of nuclear matter can be calculated. One obtains the result: 

= 5p^|«o[(l + 4)(l-2ta2) + 10x„] 

^2 

+^toti [4(1 +a;oa;i)(51 -221n2) + 7(a;o + a;i)(57-41n2)] 

H— -t\ {1 + xl)( _ 46 In 2 ) 4- xi (607 — 48 In 2) 

+^tl [ 5(1 + xl) - 78 In 2 ) + a;2(3187 - 624 In 2) 
k‘j 

+ yrOO(2 - X 0 + X 2 - 2xoX2){4:\n2 - 7) 

6 

2k^ 

+ -^tit2i2 - xi+ X2- 2xiX2)il2\n2 - 31) 

63 

10)^4 'j 

+^M''o(227-961ii2)|, (19) 

which also includes sp-interference terms proportional to toO and tit 2 - When treating the spin-orbit 
interaction at second order, the spin-trace over three polarized medium-insertions eq.(7) singles 
out the third component of the vector b = qi x q 2 , but by spherical symmetry of the integration 
region one can make replacement 63 —)■ &^/3. A good check of the intricate second-order calculation 
of asymmetry energies is provided by the relation: 

Mk^ 

- ^H«(43 + 24 In 2), (20) 
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whose first part follows by working with the equivalent interaction formulated in terms of the 
isospin-exchange operator P^- = {1 + fi-f 2)/2 instead of Pa = (1 + <?i-?2)/2. Obviously, the 
spin-orbit interaction produces different responses to asymmetries in isospin-space and spin-space. 

Finally, by reweighting the Hartree and Fock diagrams (with factors 1/2 and 1) one obtains 
the following second-order contribution to the spin asymmetry energy of pure neutron matter: 

= ^|-iS(l-x„)"(2 + ln2)-^(„(,(l-x„)(l-ii)(13 + 41ti2) 

- 3:i)"( 130 + 331n2) + ^i=(l + 12^(359 - 781n2) 

+^^ 0 ^ 2(1 — Xo)(l + X2)(41n2 — 7) + ^^l^2(l—2^i)(l + 2:2)(121n2 — 31) 

( 21 ) 

This expression includes sp-mixing terms and it depends on the expected parameter combinations. 


10^ 

189 


fFo^(227- 961n2) 


4 Effective nucleon mass and Fermi-liquid parameters 


In this section, we study the effective nucleon mass and the isotropic Fermi-liquid parameters. 
Through interactions with the dense nuclear medium, nucleons acquire a momentum- and density- 
dependent self-energy and thereby turn into quasi-particles. Long-lived quasi-particle excitations 
exist only in the vicinity the Fermi surface \p\ = kf and they are (kinematically) characterized 
by an effective mass. The ratio between the free nucleon mass M and the density-dependent 
in-medium nucleon mass M*{kf) can be written in the form: 


M 

M*{kf) 


2M2 


+ Rikf) , 


( 22 ) 


where the second term —k‘j/2M^ stems from the relativistic correction —p‘^/8M^ to the kinetic 
energy. The effects of interactions are represented by the (dimensionless) quantity R{kf) = 
—2Mkffi{kf)/3Ti‘^ with fi{kf) the spin- and isospin-independent p-wave Landau parameter [8]. 
The hrst order contribution from the Skyrme interaction Vsk to the quantity R{kf) reads: 




Mk] 

127r2 


\ 3 ti ^2(5 + 4x2)] . 


(23) 


In order to determine the corresponding second-order contribution, one has to evaluate the hve 
direct and crossed one-loop diagrams (and their mirror-reflected partners) for forward scattering 
of two nucleons on the Fermi surface |pi| = |p 2 | = fc/, shown in Fig. 2. The p-wave projection 
is performed by angular averaging with a weighting factor 3pi • p 2 • After a somewhat lengthy 
calculation one arrives at the following result: 




M‘^k‘^ f k‘^ 

+ ^o)(l - 7In2) -h -f 0 : 02 : 1 )(47 - 104In2) 

+ 3 ^i?(l + x?)(911 - 9421112) + + 8 x 2 + 5xi)(313 - 4261ii2) 


(24) 
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Figure 2: Direct and exchange diagrams representing the quasi-particle interaction at second order. 
Exchange diagrams need to be multiplied by the negative exchange operator —Pt-P„. 


which interestingly depends on the same combinations of parameters + Xg), toD(l + 

tl{l + xf), ^ 2(5 + 8 x 2 -l- 5 x 2 ) and Wq as the second-order energy per particle written in 

eq.(16). 

The isotropic (i.e. angle-averaged) quasi-particle interaction Po{kf) on the Fermi surface 
IPigl = ^/ can be decomposed into four density-dependent (s-wave) Fermi-liquid parameters: 


J^oikf) = fo{kf) + fo{kf)Ti-T2 +9oikf)ara2 + gQ{kf)Ti ■T2^r ^2- (25) 


The importance of the isotropic Fermi-liquid parameters fo{kf), f'^^kf) and goikf) stems from the 
fact that these obey exact relations to specific nuclear matter properties. The Landau relation [S] 
to the nuclear matter incompressibility K{kf) reads: 


K{kf) = kj 


d^Ejkf) 

dkj 


+ Akf 


dEjkf) 

dkf 


M 


2M2 


+ Eikf) -F 


2Mkf 



(26) 


where the expression on the right hand side is built up in the order: free nucleons with relativistic 
corrections, non-interacting quasi-nucleons, and finally interacting quasi-nucleons. The spin- and 
isospin-independent Fermi-liquid parameter fo{kf) receives a first-order contribution of the form: 


“ FT + -^^2(5-^ 4x2)] , 


(27) 


and the second-order contribution derived from the one-loop scattering diagrams in Fig. 2 reads: 

/o(fc/)^^“^^ = ^^^tl{l + xl){2 + \n2) +^toti{l + xoXi){Al + 8\n2) 

^4 ^4 

-|-—^t^(l -|- x^)(32 -|- 3 ln2) -|- 77 ^t 2(5 -|- 8 x 2 -|- 5 x 2)(8 -|- 3 ln2) 

15 63 

+^IV|(463 + 241,12)}. (28) 


Note that in /o(/c/)*'^“‘^^ the same combinations of Skyrme parameters appear as in eqs.(16,24). The 
validity of the Landau relation eq.(26) between E{kf), R{kf) and fo{kf) is immediately verihed 
at first and second order and thus serves as an important consistency check of our calculation. 
There exists a Landau relation [S] for the isospin asymmetry energy A{kf) of nuclear matter: 


Mkf) 


6M 


1 - 


^ + 


2Mkf 



(29) 
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which reflects again the composition; free nucleons with relativistic corrections, non-interacting 
quasi-nucleons, and interacting quasi-nucleons. The isospin-dependent Fermi-liquid parameter 
foikf) receives a first-order contribution of the form: 


+ 2a:o) + “^[^ 2(1 + 2 x 2 ) - ti(l -h 2a;i)] , 


(30) 


and the second-order contribution derived from the one-loop scattering diagrams in Fig. 2 reads: 


foikf) 


(2nd) _ 


Mkf J ^2 

kj 


to[(l + xl) In 2 - 2xo 


+ a;oa;i)(9 -h 32 In2) -F (xq + a;i)(41n2 - 57)] 


30 

.2 


315 


[(1 + x?) (65 + 111 In 2) + Xl (48 In 2 - 607)] 


k^ 

+ ^tl [(1 -t- xl){19 + 15 In 2) + 2 : 2(43 + 24 In2)] 

Do 

k‘j 

+ ^^ 0 ^ 2(2 + Xq + X 2 + 2a:oa:2)(41n2 — 7) 

2kj 

+ ^ 1 ^ 2(2 -|- Xl -f- 2:2 -|- 2 x 1 X 2 )(12 In 2 — 31) 

315 

k^ 

+Jn,’j(103 + 241n2) 


(31) 


It is important to note that the Landau relation eq.(29) between A{kf), R{kf) and foikf) is valid 
at first and second order. 

In an analogous way the spin asymmetry energy S{kf) of nuclear matter is related to spin-spin 
part of the (isotropic) quasi-particle interaction; 


S{kf) = 


k 


f 


6M 


V , TDfUA , 2^^/ 


+ Rikf) + 


9o{kf) 


(32) 


2M2 - TT 

The spin-dependent Fermi-liquid parameter go{kf) receives from Vsk the first-order contribution: 


y ~ 1) + ^ [fi(2xi - 1) -h t 2 (l 2 x 2 )] — , 


(33) 


and the second-order contribution derived from the one-loop scattering diagrams in Fig. 2 reads: 


Mkl 


+ 241n2), 


(34) 


with /o(A:/)*'^“^^ written in eq.(31). The Landau relation eq.(32) between S{kf), R{kf) and go{kf) 
is just as well satisfied at first and second order. 

Finally, there is the isovector spin-spin Fermi-liquid parameter g^lkf). The corresponding 
contributions from the Skyrme interaction Vsk at first and second order read: 


»;(*:/)'“> = -J + f to - il), 


(35) 
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9oikfY 


= ^^|to[31n2 - a;o(2 + ln2)] + + 321n2 - a;oa;i(41 + 81n2)] 


2 

A 

10 




f .2 


105 


tl [65 + 111 In 2 - 7a;?(32 + 3 In 2)] 


^4 ^2 

+ ^^2 [l9 + 22^2 + 8^2 + 3(5 + 8^2 + X 2 ) In 2j H—~^o^2(l + 2^2)(4 In2 — 7) 

4^4 ^4 

+ ^M2(1 + a;2)(121n2 - 31) + ^W^{103 + 241n2) 

1UO 


( 36 ) 


Note that sp-mixing terms are present in the three spin- or isospin-dependent Fermi-liquid param¬ 
eters. 

Additional checks of our calculation are provided by the Fermi-liquid parameters /no(^n) and 
Onoikn) of neutrons in pure neutron matter. In order to obtain these from the computation of the 
isotropic quasi-particle interaction Jq)(^/) in isospin-symmetric nuclear matter, one merely has to 
set Pr = 1, replace kf —)■ kn, and multiply the left diagram in Fig. 2 (with a closed nucleon ring) 
by a factor 1/2. The p-wave Landau parameter fni{kn) enters the effective neutron mass ratio 
Rnikn) with a factor —Mkn/3n'^. For the sake of completeness we list the relevant expressions for 
the neutron quasi-particle properties at first order: 

Mk^ 

RniknY^^^^ = ^ [tl(l - Xi) + 3t2(l + X 2 )] , (37) 

/no(^n)*'^^*^ = —(1 ~ Xq) + [tl(l — Xi) + 3t2(l + ^^ 2 )] ; (38) 

5'no(^n)*'^^*^ = ^(^0 ~ 4“ ^ — 1) -|- ^ 2(1 + 3 ^ 2 )] , (39) 

and at second order: 

/Vf2p2 r p2 

-a)o)^(l - 71n2)-h ^to^i(l-a;o)(l-a;i)(47- 1041n2) 

+ ^tl{l-xif{<dll-<dA2\n2) + ^tl{l+X2f{3l3-A2Q\n2) 

2k^ 

+ Y^IFo^( 881 - 1032 In 2) 



fnoiknY^^^^ 


gnoiknf^'^^ = 


Mk f k‘^ 

= - xo)\2 + ln2) + ^toti(l - a;o)(l - a;i)(41 + 81n2) 

+ ^tjil-xif{32 + 3\n2) + ^tlil + X2fi8 + 3\n2) 

15 7 

2p4 ') 

+ Y5|^o(463 + 241n2)|. (41) 

Mk f 

- xo)\\n2- 1) + ^toti(l - a;o)(l - a;i)(31n2 - 4) 

p4 op4 

+ - XiY{30 ln2-53) + ^ 3.^)2(3 ^ 2 In 2) 

4"Tn^0^2(1 — a;o)(l + a;2)(41n2 — 7) -|- ytit2(l — a;i)(l -t- a;2)(12 ln2 — 31) 

10 105 

2p4 'I 

+ ^IF2(463 + 241n2) L (42) 


11 



The pertinent Landau relations to the incompressibility Kn{kn) and spin asymmetry energy Sn{kn) 
of pure neutron matter (analogous to eqs.(26,32) but with a prefactor Mkn/'n''^ at the last sum¬ 
mand) are satished at hrst and second order. The quasi-particle interaction of neutrons in pure 
neutron matter, taking into account also the non-central tensor and cross-vector terms, has been 
calculated in ref. [U] on the basis of chiral two- and three-nucleon potentials in many-body per¬ 
turbation theory up to second order. 

5 Tensor interactions 

The Skyrme interaction Vsk in eq.(l) depends on seven parameters. This are two less than the nine 
low-energy constants Cs^t and Cj ,j = l ,..., 7 appearing in general NN contact-potential of chiral 
effective held theory [ini[iii[i2] up to order p^. A complete matching a achieved by supplementing 
to the Skyrme interaction two additional tensor terms of the Galilei- and time-reversal-invariant 
form: 


Vten 


U [3 (Jl • (fout ^2 • <fout + 3 (Ti • (^n <?2 ■ gin - <?i ■ ^2(ft^ut + 

TG [3 (Ti • ^ 2 ‘gin + 3 cfi-gin( 72 ' Qout — 2 ai ■ ^2 ^ut • ^n] • 


(43) 


Note that for both terms in V(en the spin-matrices cr{al are contracted with symmetric and traceless 
tensors depending on the external momentum-differences ^n and gout. Because of these properties 
the hrst-order contribution from Iden as well as all its interference terms with Vsk vanish, regarding 
energies per particle and the central quasi-particle interaction in nuclear matter. 

Using dimensional regularization (as discussed in the beginning of section 3) and evaluating the 
second-order Hartree and Fock diagrams in Fig. 1 with the tensorial contact-interaction Vten, one 
obtains the following contributions to the energy per particle of isospin-symmetric nuclear matter: 




Mk^f r t2 'I 

^5^1^(1525 - 1291ti2)+ «2(19-31n2) j , 


(44) 


and to the energy per particle of pure neutron matter: 

B„(i„)'“> = ^i?(19-31ii2). (45) 

The asymmetric medium-insertions introduced in eqs. (4,7) lead to the following second-order 
contributions from Uen to the isospin asymmetry energy of nuclear matter: 


. , Mk^f r „ „/ 317\ 

21(fc,.)(‘“) = t|(37l9 - 7081n2) - i^p001n2 + — j + 24* 445(121112 - 31) 


(46) 


and to the spin asymmetry energy of nuclear matter: 

Mk^ ("114-2 j-2 'I 

A(A:/)(*“) = ^(13-48 In2)+ ^(395 - 1104In2) L (47) 

Interestingly, the f 4 f 5 -interference term is absent in and The actual calculation 

of the second-order Hartree diagram produces in both cases an integrand H = Gfsdgil'^ — |g 2 d) 
for the triple Fermi-sphere integral in eq.(15), but due to the symmetry qi y-)- g 2 of the other 
ingredients it integrates to zero. By reweighting the Hartree diagram with a factor 1/2, the 
second-order contribution from Uen to the spin asymmetry energy of pure neutron matter follows 
as: 

Mk^ 

= ^^^ 5(395 - 11041n2). (48) 
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Vsk 



Figure 3: Iteration of the zero-range Skyrme interaction with one-photon exchange. 


One observes that the neutron matter quantities written in eqs.(39,42) depend only on the param¬ 
eter This feature is explained by the fact the first tensor term in eq.(37) proportional to is 
responsible for the mixing between the ^S'l-state and ^Di-state with total isospin 7 = 0. 

We continue with presenting the second-order contributions from Iden to the effective nucleon 
mass ratio: 




' 1/ f 


t^p(929 - 


1267r4 1 3 


In2) -F ^(5147 


• UOU^ Hi ^) 


(49) 


and to the isotropic Fermi-liquid parameters: 


. , Mkh f .,/2195 \ o 1 

1^4 + 88 In 2 ) + tl(2029 + 312 In2) |, (50) 

Mk^ f 32 1 

3^2(40 In 2 - 73) + ^^ 4 ^ 5(12 In2 - 31) + ^^(1081 + 408 In2) |, (51) 

Mk^ (" 4-2 'j 

9o(*:/)'“’ = -^{4(881 + 9841112) 4- (?(19 4- 72ln2)|, (52) 

Mk^ r 32 1 

^'(A)^)h-) = t2(4oin2 - 73) + -^^ 4 ^ 5 ( 8 ! - 121n2) - -|(1613 + 24241n2)|. (53) 

Note that t 4 f 5 -interference terms are present only in the isovectorial quantities fo{kf) and g'oikf). 
The simultaneous validity of the three Landau relations in eqs.(26,29,32) provides an excellent 
check of our second-order calculation with the zero-range tensor interaction Iden- The second- 
order contributions from Iden to the neutron quasi-particle properties are obviously proportional 
to tl and they read: 

M^k^ 

Rniknf^'^^ = ^^t^(5147- 65041n2), (54) 

Mk^ Mk^ 

= J^tl{2029 + 3121n2), + 72hi2). (55) 


6 Iteration of Vsk with one-photon exchange 

The Skyrme interaction Vsk in eq.(l) can be viewed as a suitable phenomenological parameter¬ 
ization of the strong nuclear force. In nuclear matter (or hnite nuclei) it affects protons and 
neutrons in the same way. This isospin-degeneracy gets lifted by taking into account the Coulomb 
interaction between the protons. A further charge-symmetry breaking contribution is provided by 
iterating the (strong) Skyrme interaction Vsk with the (electromagnetic) one-photon exchange, as 
depicted by the three-loop diagram in Fig. 3. In the case at hand the iterated diagram with two 
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medium-insertions does not vanish in dimensional regularization due to the presence of a photon- 
propagator [m^ -|- (/ ± with an inhnitesimal photon mass. The limit —)■ 0 has to 

be taken after performing the d^Z-loop integral, which results in: arctan(2|g I/m..),) —)■ 7 r/ 2 . Af¬ 
ter combination with the Pauli-bloeking corrections (supplied by the iterated diagram with three 
medium-insertions) one obtains the following contribution to the energy per proton: 


Ep(^kp) 


to(l — a^o) -2 In 2^ H- f- ti(l — Xi)(— -3 In 2^ 


f 

lOvr^ \ 

-hf2(l + X 2 ) - TT^ - 51n2j -F 


2ju. 


3M2 


14 [ ' ^V 4 

to(l-a;o)(ll-21n2) 


(56) 


where a = e^/ 47 r = 1/137.036 denotes the electromagnetic hne-structure constant and the proton 
Fermi momentum kp is related to the proton density by Pp = fc^/dvr^. Remarkably, the result 
for Ep{kp) is infrared-hnite, since an infrared regulator term lrL{kp/m-y) drops out in the end by 
angular integration: J^^dy [/(y)— /(—?/)] = 0. The last term in eq.(56) comes from iterating the 
s-wave contact part to(l — xq) with the magnetic interaction ~ where pp = 2.79 denotes the 

proton magnetic moment in units of the nuclear magneton e/2M. For the sake of completeness, 
we give also the contribution from the iterated diagram in Fig. 3 to the spin asymmetry energy of 
protons: 


Sp{kp) = ^^|^|to(l-a:o)(^ + l-21n2)+ ^ ti(l-a;i)(^^-^ + ln2) 

—^ 2(1 + X 2 ) 5 ln2j -^to{xo — 1)(2 In 2) 


(57) 


As an aside we note that the one-photon exchange with the magnetic coupling gives also (small) 
hrst-order Fock contributions of the form: Ep{kp) = —3Sp{kp) = a p^k'^/{QnM‘^). 

The charge-symmetry breaking contribution ppEp{kp) to the energy density of protons in nuclei 
could play a role in resolving the Nolen-Schiffer anomaly [91 Eg. This anomaly refers to the obser¬ 
vation that the binding energy differences of mirror nuclei cannot be explained by the Coulomb 
interaction of protons alone, using realistic nuclear wave functions. In a systematic study by Brown 
et ah na it was found that the anomaly can be resolved by either dropping the attractive Coulomb 
Fock-term Ep{kp) = —3akp/A7T or by introducing a charge-symmetry breaking delta-force. The 
hrst term in eq.(56) proportional to k^ has been investigated in ref. [18] with respect to its ca¬ 
pability of compensating the attractive Coulomb Fock-term. Such a study should be extended 
to the iteration of the complete Skyrme interaction Vs^ with the one-photon exchange. Actually, 
the novel charge-symmetry breaking contribution Ep{kp) written in eq.(56) should be included 
routinely in Skyrme calculations of medium-mass and heavy nuclei. 


7 Provisional parameter fits 

In this section we perform some provisional parameter hts in order to learn about possible improve¬ 
ments through the inclusion of second-order contributions from Vsk- We restrict the discussion to 
the nuclear matter saturation curve in the density region 0 < p < 2po = 0.32 fm“^. The equation of 
state provided by the phenomenological SLy5-interaction [19] is taken as a benchmark. Out of in¬ 
terest we perform hrst a ht to the energy per particle at hrst order: E{kf) = E{kfY^^'^'> —a^kj+a^kj. 
The empirical nuclear matter saturation point (fcjo = 263 MeV, Eq = —16MeV) hxes the two free 
parameters to the values: 03 = 30.81 MeVfm^, 05 = 8.356 MeVfm®, and the resulting saturation 
curve is shown by the full line in Fig. 4. It is considerably stiher than the realistic equation of 
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Figure 4: Saturation curve of isospin-symmetric nuclear matter. The dashed line stems from the 
SLy5-interaction [12]. The full curve is a first-order fit without the density-dependent term 

state (represented by the dashed line in Fig. 4) as evidenced by the high value K = 307 MeV of 
the nuclear matter incompressibility. While the shape of the full line in Fig. 4 can be considered 
as acceptable, there is serious problem behind it, namely the unrealistically low value of effective 
nucleon mass M*{kfo) = 0.39M, which is completely determined by the parameter a^. This (un¬ 
wanted) rigid correlation at first order is usually avoided by introducing the density-dependent 
term 1 ^ 3(1 + 

Next, we include in the fit the complete second-order contribution written 

in eqs. (16,44). The nine parameters are determined by minimizing the squared deviations from 
the realistic nuclear matter saturation curve at about 40 equally-spaced points in the density 
region 0 < p < 2po = 0.32 fm“^. The resulting fit values are: to = —3.535 fm^, ti = 1.156 fm^, 
t 2 = —3.942 fm"^, and the others very close to zero: a:o,i ,2 = 0, ^ 4,5 = 0, Wq = 0. An inspection of 
the full line in Fig. 5 reveals that a good reproduction of the realistic nuclear matter saturation 
curve cannot be obtained. The reason for this failure lies presumably in the too strong density- 
dependence kj ~ of several second-order contributions. But interestingly with the above 
parameter values, the isospin asymmetry energy A(fc/o) = *S'(fc/o) = 33.0 MeV comes out close to 
its empirical value. The previous free fit gave zero even for the spin-orbit parameter Wq. If one 
fixes it to its (large) empirical value Wq = 126MeVfm^ = 0.638 fm^ [19], the other parameters are 
fitted to: to = —3.849 fm^, E = 1.778 fm"^, ^2 = —3.341 fm"^, and a:o,i ,2 = 0, ^ 4,5 = 0. The resulting 
equation of state is shown by the dashed-dotted line in Fig. 5 and one observes that its shape 
remains almost unchanged in comparison to the free fit. The values for the isospin asymmetry 
energy A{kfo) = 36.0 MeV and the spin asymmetry energy S'(fc/o) = 32.3 MeV have changed a bit. 

From these findings one can conclude that the inclusion of the (ad hoc) density-dependent 
term |t 3 (l + is indispensable for a realistic description of nuclear matter in the Skyrme 

approach. Formally, this extra term is taken into account by making the parameter substitutions: 
to —t to + |t 3 p^t® and to^o —)■ toXo + Ito^op^t®, in all first- and second-order expressions. We have 
convinced ourselves that after this extension a good reproduction of the realistic nuclear matter 
saturation curve is possible at second order. This is in agreement with the findings of refs. HE], 
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where however incorrect expressions have been used for the second-order contributions. Certainly, 
the htting of parameters of the Skyrme interaction Vsk, when treated to second order, should go 
beyond the mere saturation curve and include as constraints all possible available properties of 
nuclear matter. Such a comprehensive analysis goes beyond the scope of the present paper. But 
the analytical calculations of a large variety of nuclear matter properties to hrst and second order 
as performed in this work lay the foundation for such a future study. 



Figure 5: Fits to the nuclear matter saturation curve with inclusion of second-order contributions 
from the Skyrme interaction Vsk- The (ad hoc) density-dependent term is left out. 


Appendix A: In-medium loop featuring s-wave and p-wave 
mixing 


<KAAA/>» 

<WW\i> 


A n 


n 


P 


P 


Figure 6: In-medium loop for an interacting proton-neutron pair. The short double-lines symbolize 
Pauli-blocking effects in the nuclear medium. 

In this appendix we analyze under which conditions a mixing between s-wave and p-wave 
components of the two-body interaction (with different parity) can occur in a fermionic medium. 
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The two diagrams in Fig. 6 represent the in-medium loop for an interacting proton-neutron pair 
in nuclear matter. The short double-line symbolizes the medium-insertion (i.e. the Pauli-blocking 
effect) for the neutron on the left side and for the proton on the right side. If the lower and upper 
contact-vertex are of s-wave type and p-wave type, respectively, the real part of the in-medium 
loop is given by a principal-value integral of the form; 

fwr - F+- IP+ri)} 

= + (58) 

where kn,p = kf{l ± 5)^/^ denote the Fermi momenta of neutrons and protons, and j = 0,1 are 
exponents of /^. The expression in the second line follows by substituting /—)■—/ in the second 
summand. It makes evident that sp-mixing can occur only in a medium with unequal neutron and 
proton densities, kn ^ kp. For small isospin asymmetries 6 the available phase-space becomes a 
(thin) spherical shell of thickness 2kfS/3. The net result for the in-medium loop with sp-mixing, 
expanded to linear order in 6, is given in the third line of eq.(58) and it involves the logarithmic 
functions: 

/o(s, k) = ^|4s-F(s^-Ffi;2-l)ln||^-^|jj—^|, Ii{s, k) = 8 + k^Io{s, k) , (59) 

where s = \P\/kf and k = \q\/kf are two dimensionless variables. In fact, the same asymmetric 
configuration as considered above is present in the calculation of the isospin (or spin) asymmetry 
energy from the second-order Hartree and Fock diagrams. This explains why sp-interference terms 
proportional to to ^2 and tit 2 do occur in the expressions A{kfY‘^^'^\ and 

written in eqs.(18,19,21). In the case of the energy per particle, E{kf) or En{kn), one has from the 
outset just one single Fermi momentum, kf or /c„, and therefore sp-mixing terms cannot emerge 
for these quantities at second order in many-body perturbation theory. 


Appendix B: Interaction through long-range potentials 


In this appendix we derive the first-order Fock and second-order Hartree and Fock contributions 
as they arise from several long-range interaction potentials. The energy per particle and asym¬ 
metry energies are calculated for many-fermion systems with either four or two internal degrees 
of freedom (e.g. isospin-symmetric nuclear matter or pure neutron matter). Our first candidate 
is an interaction potential in coordinate-space with an inverse square and an inverse fourth-power 
dependence on the particle-distance: 


hlong 


^2 /^4 

In-fgp In-fg 


(60) 


where the parameter (32 is a length and (3^ is a volume. The calculation of interaction contributions 
in the fermionic medium requires the Fourier-transform of this long-range potential: 


hlong 


2tt'^I32 

|<?out Qinl 


/^41 •?out Qin I ) 


(61) 
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where ^ut —is the Galilei-invariant momentum transfer. Taking Gong as an interaction-vertex 
in momentum-space, the evaluation of the hrst-order Fock diagram and the second-order Hartree 
and Fock diagrams in Fig. 1 gives for the four-component system (nuclear matter): 


m) = 


l32k'j iS^Mkj 


+ 


5 35 

^2^4Mkj 


+ 


2 • TT^ 

,2G - 

5 V 3 


21n2 


iSjMk^r / 731 

35 + 


81n2j , 


(62) 


A{kf) = S{kf) = 


132k j 


+ 


2/34kj 


9 15 

A^2PAMkj 


+ 


f3iMk 


9 


^^2G+^ 


45 


(1 -3G + 21n2) + 


6 

(3lMk) 

105 


21n2 


(5G + 


103 


81n2j, (63) 


where G = fi^dxx~^ arctanx = ~ 0.9159656 denotes the Catalan number. 

Note that G always results from the d^Z-loop integral belonging to the second-order Fock diagram 
with two medium-insertions. The analogous calculation in the two-component system (pure neu¬ 
tron matter) gives the following contributions to the energy per particle and the spin asymmetry 
energy: 


h/ji ( kji ) 


I32k‘^ ^ 3f34k^ P^Mk^ / . _ n 


+ 


5 35 

2/32/34Mkf, 


+ 


10 


2 2 

4G-^-l-21n2 

o 


35 


(In2 - 1 - 3G) + ^|^(5G + - 41n2) , 


(64) 


Sn{kn) = 


/ 32 k‘^ 2/34/c^ jd^Mk, 


+ 


9 15 

l32l34Mk^ 

4^ 


TT 


(11 - 12G + 4 ln2) + (5G 

105 V 


4 In 2) 
36 J 


(65) 


It is worth to mention that the (singular) /52-interaction provides an infrared-hnite result at second 
order. Ultraviolet divergences have been treated in dimensional regularization, using the rule: 
J^dl = 0. Note that the contributions of second-order Hartree diagrams (with two or three 
medium insertions) are in nuclear matter a factor 2 larger than in pure neutron matter. The 
corresponding pieces carry (besides a rational number) the coefficients 7r^/3 and In 2. 

As a second example for a long-range interaction we treat the familiar Van-der-Waals potential: 


hvdw 


12H 

\r 1 -r 2 


Uvdw — TJ" B\ Q'out <?in| i 


( 66 ) 


where the factor 12 has been introduced in order to get a simpler Fourier-transform. The many- 
body calculation at hrst and second order with the cubic interaction-vertex leads 

to the following results for the four-component system: 


SBk^f B^Mkf flllbm?, 

E(kf) = -^ +^ +27G-1281n2 

^ 63 2145 V 1050 

2B/32Mkj 


315 


BBAMk^f / 3133\ 

(30G + 53- 161n2) + —- ^(l61n2-7G-) 

^ ^ 231 V 30 / 


(67) 
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A{kf) = S{kf) = 


8Bkj B^Mkf /613343 


+ 


/ 


21 297 V 1050 

AB(32Mk^f 


+ 27G- 1281n2 


-(45G - 8 - 241n2) + 


4:Bf3iMk^f 


315 ' ' 189 

and for the two-component system: 

- ^ 8Bk^ , B^Mkl^ fU311 , 27G-641n2 

En{kn) - ~ + 195 ^ “ 


161n2 - 7G 


1369\ 


) ,( 68 ) 


63 

4Bl32Mkl 

945 


(45G-4 - 121n2) + 


11 

B(3,Mkl 
231 V 


;ln2-7G 


1141\ 


30 


(69) 


Sn{kn) = 


8Bkl B^Mkl^ /6701 27G - 64 In2 


+ 


21 ' 27 

Bf32Mkl 




350 


+ 


11 


(180^ - 191 - 48 In 2) + ^8 In 2 - 7^ - —) 

315 ^ ^ 189 V 30 / 


(70) 


The expressions in the second lines of eqs. (67-70) stem from the interference terms of the Van-der- 
Waals interaction V^dw with the (more long-range) (52- and /54-interactions. Finally, one notices that 
all interaction contributions derived here are in the form of even powers of the Fermi momentum, 
kf or kn- It is expected that the analytical results compiled in this appendix could be useful for 
certain condensed matter systems. 
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